Bulk Burning Rate in Passive - Reactive Diffusion. 



Peter Constantin Alexander Kiselev Adam Oberman 

0> ' Leonid Ryzhik 

Department of Mathematics 
University of Chicago 
(N ■ Chicago IL 60637 

P^'- February 1, 2008 

< 



o 



Abstract 



We consider a passive scalar that is advected by a prescribed mean zero divergence- free velocity 
field, diffuses, and reacts according to a KPP-type nonlinear reaction. We introduce a quantity, 
the bulk burning rate, that makes both mathematical and physical sense in general situations 
^ ■ and extends the often ill-defined notion of front speed. We establish rigorous lower bounds for 

fSJ i the bulk burning rate that are linear in the amplitude of the advecting velocity for a large class 

CO ■ of flows. These "percolating" flows are characterized by the presence of tubes of streamlines 

connecting distant regions of burned and unburned material and generalize shear flows. The 
bound contains geometric information on the velocity streamlines and degenerates when these 
I oscillate on scales that are finer than the width of the laminar burning region. We give also 

^\ • examples of very different kind of flows, cellular flows with closed streamlines, and rigorously 

I prove that these can produce only sub-linear enhancement of the bulk burning rate. 

> 

1 Introduction 

Quite often mixtures of reactants interact in a burning region that has a rather complicated spatial 
^ I structure but is thin across. This reaction region moves towards the unburned reactants leaving 

■ behind the burned ones. When the reactants are carried by an ambient fluid then the burning 

rate may be enhanced. The physical reason for this observed speed-up is believed to be that fluid 
advection tends to increase the area available for reaction. 

Many important engineering applications of combustion operate in the presence of turbulent 
advection, and therefore the influence of advection on burning has been studied extensively by 
physicists, engineers and mathematicians. In the physics literature one can find a number of models 
and approaches that yield different predictions - relations between the turbulent intensity and the 



burning rate (H, 21, M, hM). These results are usually obtained using heuristic models and physical 



reasoning. For a recent review of some of the physics literature we refer to |28, 30 1. 



The key question we wish to address is: what characteristics of the ambient fluid flow are 
responsible for burning rate enhancement? The question needs first to be made precise, because the 
reaction region may be complicated and, in general, may move with an ill-defined velocity. 

In this paper we will define in an unambiguous fashion a quantity V representing the bulk 
burning rate. V makes both mathematical and physical sense in general; we study its relation to 
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the advecting velocity field in a simple model. We provide explicit estimates of V in terms of the 
magnitude of the advecting velocity and the geometry of streamlines. We are mostly interested in the 
regime where the advection is strong but our estimates are valid for all values of physical parameters, 
and do not involve any passage to limit. They are also valid for certain advection velocities without 
symmetry. In situations where traveling waves are known to exist, the estimates we derive provide 
automatically bounds for the speed of the traveling waves. 

The main result of this paper is the identification of a class of flows that are particularly effective 
in speeding up the bulk burning rate. We call these "percolating flows" because their main feature 
is the presence of tubes of streamlines connecting distant regions of burned and unburned material. 
For such flows we obtain an optimal linear enhancement bound 

V>KU 

where U represents the magnitude of the advecting velocity and K is a proportionality factor that 
depends on the geometry of streamlines but not the speed of the flow. Other flows and in particular 
cellular flows, which have closed streamlines, on the other hand, may produce a weaker enhancement. 

We will take an analytic approach. Numerical work concerning our results will be published 
elsewhere. We consider a well-established simple model, the passive-reactive diffusive scalar equation 

r)T ?;2 

— + u{x, y, t) • VT = kAT + ^/(T). (1) 

In this equation < T < 1 represents normalized temperature, u the advecting velocity, k thermal 
difFusivity, / the reaction and vq the laminar front speed. The advecting velocity is divergence- 
free and prescribed. No feedback oi T on u is allowed in this simple model: T is passive. The 
normalization is such that the reaction rate is This is chosen so that, in the absence of advection 
(n = 0) and given (0) below, there will exist reaction-diffusive laminar traveling wave fronts that 



move with speed at least fo |2S] 



The equation (|^), derived under assumptions of approximately constant density and approx- 
imately unity Lewis number (e.g. [^) is also used to model problems in biology (|[l3[]), chemistry, 
and has other applications (pMj) but certainly does not capture all the physical instabilities present 
in turbulent combustion (|l32|, M] )- 

The type of nonlinearity f(T) we consider in this paper is concave KPP: 

/gC^, /(0) = /(1) = 0, /'(0) = 1, r(x)<0. (2) 

The prototype non-linearity of the KPP type is /(T) = T(l — T), called KPP after pioneering work 
of Kolmogorov, Petrovskii and Piskunov . 

Other important types of nonlinearity f{T) are the Arrhenius-type 

/(r) = (i-r)e-^ 

and ignition type 

/(T) =Ofor (0,1), /(T) >OforTG (0,1), 0e(O,l). 

The mathematical study of equation (|l|) concentrated mainly on two issues: existence of trav- 
eling waves and asymptotic speed, and the homogenization regime k ^ 0. Traveling waves in one 



dimension with u = were studied in the classical works [23| and [14| with their global asymptotic 



analysis addressed later in |1£] for the ignition nonlinearity, in [y in higher dimensions, and in |37] 



for variable diffusivity and ignition nonlinearity. Traveling waves with u 7^ were shown to exist for 
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shear flows for KPP as well as for a more general class of nonlinearities 0, y, 0l- Their stability was 
established in 25|. Finally, traveling waves for periodic flows u{x,y) and ignition nonlinearity 



as well as their stability were studied in [38, 39]. Probabilistic methods were applied for the analysis 



of the KPP fronts and proof of the existence of the asymptotic speed of propagation was given in ||ll 
for periodic u ^ 0. To the best of our knowledge, until now there have been no explicit estimates on 
the speed of propagation of traveling waves or asymptotic speed of propagation with n 7^ 0, except 
for the perturbative small u result of |27]. 

The homogenization regime k — > 0, when the front width goes to zero, was extensively studied 
for KPP-type nonlinearity and for advection velocity that is periodic and varies either on the integral 
scale |jl5|, 17] or on a small scale that is larger or comparable to that of the front width 24, 



|lO| , PI , [2^ . Recently a similar result was established for random statistically homogeneous ergodic 
advection velocities [34]. A thorough review of most of these results, both for traveling waves and 
homogenization techniques is given in ]^^. The result of homogenization procedures is an effective 
equation valid in the limit k ^ 0. The effective equation is typically a non-trivial Hamilton- Jacobi 



equation [^, 11 , 24 1 . Homogenization is usually very efficient when a mean field captures the essence 
of the question asked. When that is not the case important information is lost in the limit. 

For simplicity of exposition we will consider the reaction-diffusion-advection equation (|lj) in a 
two-dimensional strip 

n = {x,y: X £ (-00, 00), ye [0,H]}, 

but the methods we introduce work in any dimension and for more general classes of domains. The 
boundary conditions are either Neumann 

Ty{x,0) =Ty{x,H) =0, (3) 

or periodic in y: 

T{x,y,t)=Tix,y + H,t). (4) 
The flow u = {ui,U2) is incompressible: 

V • u = (5) 
and has zero normal component at the boundary 

u ■ n = on dQ (6) 

in the case of Neumann boundary conditions. Furthermore, we assume that the total flow through 
the strip is zero: 

H 

ui{x,y)dy = (7) 





to eliminate the drift caused by the mean flow. We assume 1>T>0,T = 1 being the stable 
(burned) state of the system, while T = is unstable (unburned) state. The equation has maximum 
principle ]^], so if the initial data is in the [0, 1] range then the solution remains in the same range 
for all times. We also assume that the solution is localized, that is 

T(x,y,t) = 1 - 0(e^^) for x < 0, T{x,y,t) = 0{e-^'') for x > 0, (8) 

[Vr[ = O(e-^l^l) for some A > 0. (9) 

If such conditions are satisfied initially then they are valid for all subsequent times (see Section §). 

We are interested in a general situation, when traveling waves solutions may not be relevant or 
not exist. 
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We introduce a natural quantity that measures the typical burning rate 



V{t) = I Tt{x,y,t 



dxdy 
H ' 



We call V{t) the (instantaneous) "bulk burning rate" and its time average 

t 

{y)t = 7 / V{s)ds, (y)oo = liminf(y)t 

t J t^oo 


simply "bulk burning rate". Because T is non-dimensional, V has units of length per time, i.e. of 
velocity. Note that when T{x, y, t) is a traveling wave front-like solution T{x — ct, y, t), then V is the 
speed of the front, V{t) = c. But V is defined for very general initial data and more general equations 
and does not require assumptions about the nature of the burning region. The word "bulk" refers 
to the fact that we take a space (or space-time) average, capturing only bulk or large scale effects. 

Our first result shows that no matter what the advecting velocity is, it cannot slow down the 
bulk burning rate below a universal lower bound, of the same order of magnitude as the laminar 
front speed. Let us denote 

a = - inf f'{T) > 0, P = - sup /"(T) > 0. 

0<T<1 0<T<1 

(Note Q = l,/3 = 2 for KPP). 

Theorem 1 There exists a constant C > such that for any advecting velocity u{x,y,t) satisfying 
(j^ and (Q), any solution of the passive-reactive diffusion equation (jjj with boundary conditions 
^ and ^) or the bulk burning rate V{t) obeys the lower bound 




V{t) > Cvo\ - e-°^o*/2«). (10) 



One of the applications of this theorem is in a homogenization regime, where the reaction is very 
weak (see Appendix^). 

As far as the general upper bound is concerned, it is easy to show (see Section ^) for a very 
general class of velocities u{x,y,t) that if the initial data To(x,y) satisfies (^ with A = vo/2k, then 

< ^+ ll^^llloo + ^^O, (11) 

with the constant length Lq depending on the initial data Tq only. Here ||ui||oo denotes, as usual, the 
supremum of \ui\ over the whole domain. Therefore, the bulk burning rate may not exceed a linear 
bound in the amplitude of the advecting velocity. For a large class of flows we prove lower bounds 
on the bulk burning rate that are linear in the magnitude of advection. For instance, a corollary to 
Theorem 4 of Section 4 concerning mean zero shear flow of the form 

rH 

u{x,y) = {u{y),0), / u{y)dy = 
Jo 

can be stated simply as 
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Theorem 2 There exists a constant C > that depends only on the nonlinearity f but not upon 
u{y) and Tq such that, for any solution T{x, y, t) of the passive-reactive diffusion equation ^ with 
boundary conditions or and any t > tq = max[^, ^] the bulk burning rate obeys 



"0 



2 \ II Il2 

K \ \\u\\ 



(y)r>c{i + :^] ^ (12) 



0"'u , 



H 

where \\u\\i = / \u{'y)\— and hu ^ 



H ^ u' 







Recall that the normalization of the reaction rate in {|T|) is chosen so that the laminar traveling 
wave front speed is vq no matter what k is. If we allow k to vary while keeping the coefficient M in 
front of the reaction term fixed, we find that the bound is still independent of k, however the time 
To in which it is reached behaves as 

Actually, we prove a far more general geometric estimate (Section ^), from which (|l^) follows. 
This general estimate provides a non-trivial lower bound also for the case when the ratio of and 



norms of u{y) becomes small. An important feature of the bound (12) is the presence of the ratio 
of the characteristic scale hu of variations of the advection velocity and the reaction scale I = k/vq. 
The estimate degenerates in the case hu «C I which is expected from physical considerations since 
additional wrinkling on the scales smaller than the reaction scale should not accelerate the front. 

Furthermore, we consider time dependent shear flows, and obtain a similar lower bound on the 
bulk burning rate: 

u 

(V) > i^(/in,n,^)||M||oo. 



U 



Here r* is a typical time scale of the flow, defined similarly (but not identically) to hu- The prefactor 
K becomes smaller when either hu becomes smaller than the reaction scale / = k/vq, or the time 
scale T^, is faster than the reaction time Tc = k/vq or the time th = H/vq it takes the reaction to 
traverse the cross-section. The precise formulation of our result for the time dependent shear flows 
is given in Section |5[ 

Finally, in Section |6| we consider a generalization of the time independent shear flows. Namely, we 
consider "percolating flows" , flows that have two or more (sufficiently regular) tubes of streamlines 
connecting x = — oo and x = +oo. These flows are not necessarily spatially periodic and can have 
completely arbitrary features outside the tubes of streamlines. We show that the bulk burning rate 
is still linear in the magnitude of the advecting velocity, no matter what kind of behavior (closed 
streamlines, areas of still fluid, etc.) the flow has outside the tubes. The proportionality coefficient 
depends on the geometry of the flow. Thus, we identify a broad class of the flows which increase 
the bulk burning rate linearly with the amplitude of the flow, the fastest possible rate of increase. 
We also show that in general the dependence of the bulk burning rate V{t) on the magnitude of the 
advecting velocity may be sub-linear. An extreme example is provided by shear flow perpendicular 
to the front (in periodic boundary conditions) where there is no significant enhancement: the bulk 
burning rate remains uniformly bounded as the magnitude of the advecting velocity tends to infinity. 
Other examples are certain cellular flows (flows with closed streamlines): for every a £ (0,1) we 
construct cellular flows for which V{t) is bounded above by C||ii||J^ (Section 0). A comparison of 
the results of this paper with extensive numerical studies will be presented in a companion paper 
with Fausto Cattaneo, Andrea Malagoli and Natalia Vladimirova 

For the rest of this paper, we assume that the reaction f(T) satisfies (^, the initial data satisfy 
(^), and the advection velocity satisfies ||iii||oo < oo, ||Vti|[oo < oo. As we will see in the next 
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section, such assumptions on u ensure that the locaUzation @, @ of the initial data is preserved 
during the time evolution. Throughout the paper we denote by C various (not necessarily equal) 
constants which may depend only on reaction f{T). 



2 Preliminaries and an upper bound on the bulk burning rate 

Our considerations in this section follow the general ideas of Q. We show in this section that the 
boundary conditions (|8|) are conserved by evolution and establish the simple upper bound (11). 

Lemma 1 Assume that the initial data TQ{x,y) satisfies the following bounds: 

T(x,y)<Coe-^^ 1 - T(x, y) < Coe^^ |VT| < ^e-^l"l, Co > 0, A > 0. (13) 

II 

Let ci > ||ui||oo + + and C2 > \\ui\\oo + '^-^ + ^ + Then 

T{t,x,y) < Coe-^o(^-^i*), 1-T{x,y) < C7oe^»(^+^i*\ (14) 
lVr| < ^e^^(^^"2*) for any t. (15) 

Proof. The proof is an application of the maximum principle. Note that T satisfies an inequality 

Tt + u-VT- kAT - -P-T < 0. 

Introduce </>(x,t) = q-K^-^i^)^ then 

2 2 

0i + u • V0 - /tA</> - ^(j) = A(c -ui- kX- — )(/> > 
Ak 4kA 

by our assumptions on c. Applying the maximum principle to the function w = CqcJ) — T we obtain 
the first estimate in ([l^ ) . To get the second estimate we note that G = 1 — T satisfies the inequality 

Gt + u-VG- kAG < 0. 

We let then tp{x, t) = 6-^0(^+^1*) and proceed as before applying maximum principle to wi = CQip — G. 
The decay of |VT| is obtained as follows. Let P = |VTp, then P satisfies the equation 

Pt + u-VP- kAP + 2(|VT^|2 + \VTy\^) = :^f'{T)P - 2{T^u^ ■ VT + TyUy ■ VT). 

Therefore if we let K = ^ + 4||V'u||oo, then we get 

Pt + u-VP- kAP -KP <0 



and then ( |15D follows as before. □ 

Lemma |l] implies that the bulk burning rate cannot be larger than ci = vq + ||mi||oo provided 
that the initial data decays fast enough. 

Theorem 3 Assume that To{x,y) < Coe~'^^ and 1 — To(x,y) < Coe'^^ with A = — , then 

{V)t<^ + vo + \\uiU. (16) 

Vot 
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Proof. We have 

{Vm = 1 J^ds J ^T,{s,x,y) = i J ^[nx,y,t) - To{x,y)]. 

Lemma |l| implies that T{x, y, t) satisfies the bound 

r(x,y,t)<Coe-^(^-^i*), 

with ci = ||tii||oo + kA + j^. Then we have 

{V){t) < -j^ ^J^J(i-To)-{l-T))dx + -J^ -|y^ [T{x,y,t)-To{x,y)] 

+ - / T7 / Tdx<——+ci<— h||ni||oo+wo 

t Jo H Jc^t tX tvo 

with our choice of A = f^-Q 

A simple corollary of Theorem ^ is that a shear flow in the direction perpendicular to the front 
propagation does not enhance the bulk burning rate. 

3 A universal lower bound for bulk burning rate 

In this section we prove Theorem Let us integrate (||) over the set and obtain 

Integration by parts is justified for any t by Lemma |l]. A straightforward computation that uses ([l|), 
the boundary conditions (||) or (Q), (P), and the incompressibility of u{x,y) shows that 

^>— y |Vr| y. (18) 

n 

The two simple equations (|l7| ) and (jl^) are the basis of our technique for deriving lower bounds 
on the burning velocity. The temperature drops from 1 on the left to on the right. The reaction 
/(T) is large where T takes values in some region strictly between and 1. If T varies slowly, we 
get a good lower bound for V{t) from (|l^). On the other hand, if T varies fast, then the |VTp 
term in (^) is large, which will give a lower bound on the bulk burning rate. The equations ( p^ 
and ( |l8|) are thus complementary in this respect. The following Lemma gives precise meaning to 
the statement that the reactive term in ( p!7| ) and the gradient term in (^) cannot be simultaneously 
small: 

Lemma 2 Let f{T) he a function of concave KPP type ^ and assume that the continuously dif- 
ferentiahle function T{x, y) satisfies the following assumptions: 

(i) o<r(x,2/)<i, 

(a) lim T(x,y) = 1, lim T(x,y) = for every y £ [0, 1], 
Then there exists a constant C > ( independent of the function T ) such that 

j f{T)dxdy J \VT\^dxdy > CH^. (19) 
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Proof. We may assume that / \VT\'^dxdy < oo and / f(T{x,y))dxdy < oo, otherwise (|T^) is trivial. 

n n 

Let y E (0, 1) be such that 

dxdy 



J \VT{x,y)\'dx<3j |VT|2^ 



and 

oo 

dxdy 



f{T{x,y))dx<3 f{T{x,y)) 



H 



Then there exist xi, X2 such that \T{xi,y)—T{x2,y)\ > 1— e while f{T{x, y)) > Ce for all x G (xi, X2) 
because of the boundary conditions (ii) and property (0) of the non-linearity. Then we have 



Ce|xi-X2| <3 j f{T{x,y)) 



dxdy 
H 



n 

and 



n 

Multiplying these two equations we obtain 

Ceil - e?H^ 



j f{T)dxdy j \VT\^dxdy > 
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which proves Lemma □ 

Lemma (|T7|) , and (|l^) imply that 



dt Ak ~ IQkV 

Therefore 



~ 4a 



from which Theorem || follows. 

We remark that a simple variation of Lemma |2| allows to prove Theorem || for more general 
domains than a strip. 



4 Bulk burning rate in shear flows 

Let us consider the passive-reactive diffusion equation in a shear flow in two dimensions: 

Tt + u{y)T, = kAT + ^/(T). (20) 

The boundary and initial conditions are as in (P) or (Q) and (^), @. The flow u{y) is continuously 
differentiable and has mean zero (a non-zero mean can be taken into account by a simple change of 
variables): 

rH 

u{y)dy = 0. 



We prove now an estimate for the bulk burning rate which is more general than the one presented 
in Theorem 0. 
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Figure 1: The structure of the shear flow. 
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Theorem 4 Let us consider an arbitrary partition of the interval [0, H] into subintervals I 
[cj — hj,Cj + hj] on which u{y) does not change sign. Denote by the unions of intervals 

Ij where u{y) > and u{y) < respectively (see figure [ip. Then there exists a constant C > 0, 
independent of the partition, u{y), and the initial data To{x,y), so that the average burning rate 
{V)r satisfies the following estimate: 



{V)r > C 




\u{y)\'^ + c- 



E 




I ( w^y 

\u[y)\j^ 



\ 



J 



(211 



for any t > tq = max 



K H_ 




. Here I = k/vq. The constants c± are defined by 



Remarks. 1. The exact choice of the unions of intervals D± is left to us. Given a shear flow, 
we should attempt to pick in a way to maximize the bound (^). We provide a simple example 
after the proof which illustrates how this bound works. 

2. The proof of this bound becomes much easier, and also extends to more general types of reaction, 
if we assume Tt{x,y,t) > 0. This is the case for traveling waves, the existence of which has been 
proven for shear flows and various types of reaction in 0]. The minimal speed of traveling waves 
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also provides lower bounds for asymptotic speed of propagation of any front-like data We 
plan to address the results which one can prove along these lines in subsequent publications. The 
disadvantage of an a priori assumption Tt > 0, however, is that we cannot get an estimate on the 
time required to reach the lower bound from any initial data, and cannot extend the results to 
time-dependent flows, as we do in Section ||. 

Proof. The plan of the proof is as follows. We know that the bulk burning rate V{t) satisfies the 
inequalities 

V{t) + > ^ / \VT\^dxdy. (22) 



ovq dt a Jn 
and 

vit) = |i / /(D^^. (23) 

4k Jn H 

We will be able to get a bound from below in terms of u for the combination of the terms on the 
right-hand sides of these equations, which then provide a bound for {y)r- The key will be to integrate 
in X along the streamlines of the flow and then to average several times in y and t to bring Tyy and 
Tt terms to a form convenient for the estimate. 

Consider one interval of D+, Ij = [cj — hj,Cj + hj], and so that u{y) > for y e (cj — hj, Cj + hj). 
Integrating (|2^) over x G iR, we obtain 

j Ttdx-K j Tyydx-^j f{T) dx = u{y). (24) 
MM M 

Therefore 

J Ttdx — hi j Tyy dx > u{y). (25) 
M M 

(It might appear that we make the estimate rather crude by dropping the f{T) term. There is 
however heuristic and numerical evidence that in many situations this term is insignificant in the 
regions where u{y) > 0, and that the front is quite sharp in these regions. In contrast, in the regions 
where u{y) < the burning region is wider and the f{T) term will not be discarded there.) Now we 
estimate both terms on the left hand side of (p5[). Let us begin with the second derivative term. To 
reduce the order of differentiation, we employ the following averaging in y : 

hj/2 hj/2+j Cj+5 Cj+hj 

I d-f J d5 I ■dy= I Gihj,y-Cj)- dy, (26) 

hj /2-f Cj-5 Cj-hj 

where the kernel G{h, ^) can be computed explicitly as 

G(,.^) = |^('.-lfl)^-(t-l^l)^ i?i<V2 

I ^2(h-\(\r-. h/2<\(\<h 

Observe that the function G{h, has the following properties 

G{h, 0>0 for all ^ G [-h, h] 

G{h,0 < ^ for ^£[-h,h] (28) 
G{h,0>^ for 



Let us apply the averaging procedure (|26D to the equation (gq). We have 
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Lemma 3 



dx / G{hj,y - Cj)Tyy{x,y)dy 



< C 



Cj+hj 



Cj+hj 



J dx J \VT\^dy + jdx J f{T)dy\. (29) 



V K c,-h 



Proof. We are going to split x € iR into two sets. In one set, the Ly norm of the gradient will be 
large and we will use it to estimate the second derivative term. In the other set, the variation of the 
temperature will be small, and we will use the reactive term to bound the second derivative term. 
More precisely, let p be a number such that \plphj = 1/3 and define the set Vjp C Mhy 



V 



3P 



eJR : J \VT\'^{x,y)dy > phj 



so that 



dy\Ty{x,y)\<J- I \VT{x,y)\'dy = 6hj / \VT{x,y)\'dy 



for X E Pjp. Notice that for such x, according to (p6| 



hj/2 hj/2+-y 



G{hj,y-Cj)Tyy{x,y)dy= J dj J d6{Ty{x,Cj + 6) - Ty{x,Cj - 6)) 

hj/2-"i 



< 



' '"J 

3h] J \VTix,y)fdy. 



For X outside T>jp, we use the representation 

Cj+hj 

J G{hj,y - Cj)Tyy{x,y)dy = 



T{x,y)dy 



hj/2<\y-Cj\<hj 



\y-Cj\<hjj2 



We need the following crucial observation: 

Lemma 4 Assume that ^ 

\T{x,yi) - T{x,y2)\ < -, 

then we have 

\T{x, yi) - T{x, y2) \ < C{f{T{x, y{) + f{T{x, y^)) 
for any yi, y2 G {cj - hj,Cj + hj). 

Proof. Let us denote Ti = T{x,yi), T2 = T{x,y2), then we have using (Q) 

This proves Lemma |^. □ 
Notice that if x ^ "Djp we have 

Cj+hj 



\Ty{x,y)\dy < ^/2ph 



(30) 



(31) 



T{x,y)dy. (32) 
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and therefore \T{x,yi) — T{x,y2)\ < 1/3 for any yi,y2 G {cj — hj,Cj + hj). Applying Lemma § to 
(|32|), we get 

Cj+hj 



Cj+hj 

dx I G{hj,y - Cj)Tyy{x, y) dy 



<cjdx J f{T) dy. 



This completes the proof of Lemma ^. □ 

Averaging equation (p^) and applying Lemma H we obtain 



dyu{y)G{hj,y-Cj)< dx / dyG{hj ,y - Cj)Tt{x,y) 



(33) 



Cj+hj 



Cj+hj 



+c 



hjK J dx J dy\VT{x,y)\ + n J dx J dyf{T{x,y)) 



Next we consider the intervals where velocity u{y) < 0. These are analyzed similarly, except that 
now we do not discard the last term on the right side in (p^). The estimate analogous to (^) is 



Cj+hj 



dy\uiy)\G{hj,y - Cj) < ^ I dx / dyG{hj,y - Cj)f{T{x,y)) 



(34) 



Cj+hj 

dx / dyG{hj,y - Cj)Tt{x,y) 



I 



h]K^dx J dy\VT{x,y)\^ + kJ dx J dyf{T{x,y)) 



Thus we succeeded in replacing the second order derivative term with expressions directly linked 
to the burning rate. Now it remains to estimate the time derivative term. Summing ( [3^ ) and (^ 
over Ij C 15+ and Ij C respectively, and using properties (|28| ) of the kernel G, we get 



1^ J dyj dx|Vr(x,y,t)|2 + ^ J dy J dxf{T{t,x,y)) 



D+ M 



D+ M 



(35) 



and 



+ dy dx 



G{hj,y- Cj 



-1 ^j+T 



h] + P 



-Ttix,y,t)>C E 1 + 



I,CD+ 



K J dyj dx\VT{x,y,t)\^ + ^ J dy J dxf{T{t,x,y)) 



D- M 



D_ M 



\u{y)\dy 



^3 2 



(36) 



E ldyldx^^^^^y^Ux,y,t)>CY (l + i^ 



-1 ^^^^ 



\u{y)\dy 
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(recall / = h/vq). Let us choose the weights m+ and m_ according to 

G{hj,y - Cj) 



E Idy 



m±- ^ J 

Set also M = max(m+,m_). Notice that by the properties of G, we have 

1 m± 1 

— c+ < < — c+ 

16 ^ - M - 4 ^ 

for constants c± in the formulation of the Theorem. Let us define measures 



(37) 



, m±xi,{y)G{hj,y - Cj) 



MH{h] + /2 



(38) 



(39) 



Here we denote, as usual, by xs the characteristic function of the set S. Multiplying (|3^) and (pq) 
by and m_ respectively, and adding them together we obtain 

K / iVTp^ + J I /(T)^ + |dx /d.+ (y)r, -Jdx ]du^{y)T, (40) 



K 



> C 



/2 



^1 + ~ 



M 

hj 

-1 + ~ 



We have the following 
Lemma 5 For any ri, T2, 



H 



H 



J dt i J dx J dvj^{y)Tt — j dx j du_{y)Ti 



<ci:f-/ivr(.,,,r.)p^^^ 



n \K 

dxdy 



m 



+ 1 + 



f{T{x,y,Ti)) 



dxdy 



Proof. By definition (3£) of the measures v±, their total weights are equal: 



H 



H 



dv. 



< 1. 



(41) 



It is easy to construct a measure preserving bijective map ^*(y) : -D+ -D-, so that i^+(5') 
^-{^{S)). Then we can write 



T2 / H H \ 

JdtlJdxJ du+{y)Tt - j dx j du^iy)Tt (42) 



M 

H H 

j dx j dz.+ (y)(r(x,$(y),ri)-r(x,y,ri)) + j dx j dz.+ (y)(T(x, y, ra) - r(x, $(y), ra)). 

JR JR 
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Consider the first term on the left hand side of (14 2|). We spht all x E iR into two sets, saying x E S 
if there exists y such that 

\T{x,y,Ti)-T{x,^y),Ti)\ > ^ 
Using the same argument we applied in the proof of Lemma |2|, one can show that 



H 



H 



VT{x,y)\ dy f{T{x,y))dy>C supy^^y^^[o^H]\T{x,yi) - T{x,y2) 











where C is some universal constant, depending only on /. Therefore, for every x & S, 

H H ^ 



J\VT{x,y)\'dyJ f{T{x,y))\ >C, 







and hence, using the fact that the total weight of does not exceed 1 by (^T]), we have 



H 



H 



H 



du+{y)\iT{x,^y),n)-T{x,y,n))\ < C 



K 



\VT{x,y)\^dy + ^ / f{T{x,y))dy 



(43) 



For X ^ S, we have 

H 



H 



H 



J du+{y)\Tix,y,n)-T{x,^y),n)\<C J du+{y)[f {T{x,y)) + f {T{xMv)))] < C j f{T{x,y))^ 



(44) 



by Lemma Q and (|3^) . Equations ( |43| ) and ( p9| ) together imply Lemma ^. □ 

Theorem ^ now follows from Lemma ^, relations (^) and (^), and inequality (^). Given time 
interval [0,r], apply the following averaging to the both sides of (^): 



i+7 f+<5 



— J dj J d5 J dt = — I G{-,t - -)dt 

^-y 



(45) 



A direct computation using (22) and (^) shows that the left hand side of ( |40| ) after averaging does 
not exceed 

C 



(46) 



On the other hand, the right hand side of (^) is independent of time, and averaging ( ^5[ ) results in 
multiplication by a constant independent of r. Hence from (46) we see that if 

K H 

T > To = max^, — , 



then 



{V)r > C 



-1 



V 



I,CD+ 



-1 + ^ 



\u{y)\dy + 



m- 



\u{y)\dy 



r---2. 

'^3 2 
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This proves Theorem ^ (recall that m± and c± are related by (^)). Q 

Remark. At the expense of making the proof slightly more technical, the total width of the strip 
H in the formula for tq can be replaced by an often smaller value H, which is introduced as follows. 
Consider the function g(y) = i/_|_[0,y] — z/_[0,y] on [0,if]. Then H is defined as a maximal distance 
between two neighboring roots of g. It is straightforward to generalize the proof of Lemma ^ to yield 
this result, by taking a specific measure preserving function <I> which maps to only within the 



Vo 

meaning: this is the time needed to burn across the scale on which the shear flow u wrinkles the front. 



intervals between the neighboring roots of g. The characteristic time — has a clear intuitive physical 



Example. Let u{y) = uq sin . We can take intervals Ij as half-periods of u where it does 
not change sign. Factors c± are equal in this case. Then Theorem implies that for any t > tq = 



maxf-^, —1 we have 



According to the above Remark, it is easy to see that in this example, H in definition of tq can be 
replaced with H = H/n. The map <^ in this example can be taken to map half-periods where u is 
positive on the neighboring half-periods where u is negative. 

Theorem ^ is a direct corollary of Theorem 

Proof. Let us denote by IS*! the Lebesgue measure of set S. Define the sets 
•^±=|ye[0,F]: ±u{y)>\j^ |^(y)|^ = l||n||i|. 

Notice that since u is mean zero, ||u±||i = ^H^^Hi (here u± are the positive and the negative parts of 
u). Therefore 

-||u||i(i? - \T±\) + ll-ullool^il > -\\u\\iH, 

and so 

> 1;^. (47) 

Let hu = }} }\^ , then for any y € T± and any y' ^ (y — hu/8,y + hu/8) we have \u(y')\ > -||^i||Li. 

If IIloo 8 

It is easy to construct unions D± = Ujlf of non-overlapping intervals 1^ = {y^ — hu/8,y^ + hu/8) 
with y^ G J^±, such that \D±\ > — Then we have: 



, , ^,dy 1 \\u\\i , ^ , , , , 



Combining (P), (||), and (|T|) we get 




□ 
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5 Time dependent shear flows 



One may expect linear growth of the bulk burning rate in the amplitude of the flow, but the temporal 
characteristic scale t^, of the variations of the flow will play a role similar to that of the scale hu 
in Theorem |2[ That is, too rapid oscillations will diminish the enhancement of the bulk burning 
rate. Consider two systems of intervals Ij in [0, H]^ Dj^ and D^. At this point we do not make any 
assumptions regarding the behavior of u{y,t) on D±. Let us introduce the notation 

J(t,u) = c+ 5: 1 + / n(t,y)^ (49) 

'^i 2 

-1 + — ^ 




1 . 



2 



{V)to,r>C{l+[^\ \ J(to,r,u), (51) 



where c± are defined as in Theorem Given a starting time to length of the time interval r, 
we define 

J{to,T,u) = ^ J G{^,t-to-^)J{t,u)dt. (50) 

to 

We also denote (l^)to,T the average of V{t) over an interval of time of duration r starting at time to '■ 

to+T 

iy)to,r = ^ J V{t)dt. 

to 

We have 

Theorem 5 For any choice of the intervals Ij C [0, H] in D± and any r, to 
where Tn = maxl-^, —]. 

Proof. The proof is a direct corollary of the proof of Theorem ^ We remark that the choice of the 
averaging in time ( p5| ) in the last stage of the proof is not the only one possible. Given a particular 
time dependent flow at some initial moment to, one can try to adjust the averaging procedure to get 
a better lower bound. □ 

Remark. Similarly to the remark after the proof of Theorem ^, we can replace H in the definition 
of To by a smaller value H (defined in that remark). 

To clarify the meaning of Theorem ^ we make several observations and consider two examples. 
The general way one can apply this theorem is as follows. Given a moment of time to, we try to 
choose T and D± so as to maximize the lower bound (|5l|). There is a certain tradeoff involved in 
choosing r. If we take r to be small, it is likely that we can find D± so that velocity u{y, t) does 
not change sign there during time interval [to, to + t], staying positive on and negative on D_. 
Then there is no cancellation in equation ( |50|) defining J{tQ,T,u). However, the factor (1 + ^)~^ 



16 



may become very small if r <C tq. If we take r large, some cancellation is likely to occur in (50), 
making the bound weaker, unless the shear flow varies on time scales still larger than r. In the flows 
which oscillate in time on the scale smaller than tq, we will not be able to avoid either cancellation 
in (|50| ) or small factor (t/tq)^ in the bound (|5l|) , and will end up with weaker lower bound than 
we would have gotten if the flow varied slower in time. Notice that in any case the bound grows 
linearly with the amplitude of the flow, and there are factors reflecting the moderating effect of fast 
oscillations both in space and in time. 

If we want to know the average of the bulk burning rate over a long period of time, much larger 
than the typical time scale of the flow, we can use Theorem ^ by splitting this long time period into 
appropriately chosen smaller ones and getting lower bounds on the averages over these smaller time 
intervals. Combined, they will also give us an estimate on the long time average. 

Example 1. Consider a flow 

27rny 

u{y,t) = uq sm ZTTLot sm — — — . 

H 

To get an estimate on long-time average of the burning rate for such flow, consider to = 0. Set 
r = and take 

n n 
^+ = U ^2,-1, ^- = U ^2, ' 

where Ij = ( ^''^2n^ ' Then we get 



J(0,r,u) = uo(1 + 4(toc^) 



2\-l 




27rny 
sm — ^ I dy 



> C(1 + 4(to^)2)-i (^1+ ^2 

with constant C dependent only on reaction /. A similar estimate is valid for to = we only need 
to switch D±. Therefore, we get that for any to and any ti > 

{V)t,,r, > C{1 + 4(ro^)2)-i (^1 + uo. (52) 

It is not difficult to obtain estimates for averages over times smaller than ^ , but these would generally 
(and naturally) depend on the choice of starting time to- 
Example 2. Consider 

2irn(y — ct) 
u{y, t) = uo sm — . 

This is a flow which shifts in y direction. We assume for simplicity that the boundary conditions 
for T are periodic in y. Given any to, pick r = Time r is chosen so that during this time, the 
regions where u is positive and negative do not shift completely; there are regions where velocity 
stays positive or negative during [0,t]. Take 

n n 
3=1 j=l 
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with Ij = ( ; ^%nH^ )- ^ direct computation of J{tQ,T,u) shows the following bound: 

(V),,.>c(l+(?^)y' (l + !^)"„„. (53) 

where C may depend only on reaction function /. We note that in this example, it is to easy to show 
that (^3|) extends to any averaging time r, independently of the starting time to : 

{y)to,r>c(l + i] 7l + ^?l 'no. 



2 I- ' H2 



T 

6 Percolating flows 

We now consider a more general class of flows, which we call "percolating". By this we mean that 
there exist at least two tubes of streamlines of the advecting velocity u{x,y), one of which connects 
X = —oo and x = +oo, and the other one goes from x = +oo to x = — oo. More precisely, let us 
assume that there exist regions Dj' and DJ , j = I, . . . N such that each of them is bounded by the 
streamlines of u{x,y), and the projection of each streamline of u{x,y), contained in either Dj~ or 
D~ , onto the x-axis covers the whole real line (these projections need not be one-to-one, however). 
As before, we denote D± the union of all respectively. 

Our considerations in this section will follow closely the ideas of the shear flow case. However, 
there are two natural geometries in the problem. The Laplace operator is best described in Euclidean 
coordinates, while for the advection term the geometry of streamlines imposed by the flow is most 
natural. In the case of the shear flows these geometries coincide, but generally they are at odds. 
Due to this fact, additional technical difficulties arise when we consider percolating flows. 

We assume that the streamlines in are sufficiently regular, so that inside each there exists 
a one-to-one change of coordinates {x,y) — (p,0), such that p is constant on the streamlines, 
while 9 is an orthogonal coordinate for p (with a slight abuse of notation we shall use the same 
notation (p, 6) in all ^ although these coordinates are not defined globally). Moreover, u ■ V9 > 
in Dj', while u-VO < in each DJ . On D^, p varies in [c^ — h^,c^ + h^], while 9 varies in {—oo, oo). 
See figure ^ for a sketch of coordinates {p,9). The square of the length element inside each set 
is given by 

dx^ + dy^ = Ef{p,9)dp^ + E^{p,9)d9'^. 



We assume that the functions Ei^2 satisfy the following conditions. They are bounded from 
above and below: 

C-^ <Ei,2{p,9) <C (54) 



uniformly on all Dj . Moreover, the function 



ujip,9) 



E2{p, 



Eiip,9) 



satisfies the following bounds: 

< \uj{p,9)\ < C, 



C 

< -rr on Df respectively. 



(55) 



(56) 



with 2h^ being the absolute value of the difference of the values of p on the two components of the 
boundary dD^ (recall that are bounded by two streamlines of u{x,y)). 

Conditions (]5^ ) and (|5^ ) are satisfied, for instance, in the following examples: 
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Pi 

Figure 2: Curvilinear coordinates {p,0). 



± 



(1) A flow u{x,y) = Uu{x,y) with U being a scalar, and the flow u{x,y) satisfying on Dj 

<\u{x,y)\<C, \V xu{x,y)\<^, (57) 

a 

where d is the maximum length of a level set of 9 inside D^. 

(2) More generally, it is enough to ask that in each there exists a function constant on the 
streamlines of u such that 

C-^ < IWI < C, lAiPl < ^. (58) 

a 

We remark that (1) is a particular case of (2) where ip is taken to be a stream function of the flow 
u. 

We do not make any assumptions on the behavior of the streamlines of u{x, y) outside the regions 
Z?+ and In particular, there may be pockets of still fluid, streamlines may be closed, etc. (see 
Figure 1). 

Another assumption concerns the relative measure of the sets D±, which we assume to oscillate 
not too wildly. More precisely, let Rab denote the rectangle Rab = [a,b]x x [0,H]y, and let D'^ = 
D± n Rab- We define the measures 

G{hj,p- Cj) 



p±[a,b] = J dpd9Y,Ei{p,9)E2ip,e)- 



jyab 



H{hj + /2) 



with I = k/vq and the function G{h, p) defined by (pT]). Then we assume that there exists a partition 
of the real axis 

. . . < x-n < . . . < x_i < xq < xi < . . . Xn < ■ ■ ■ Xj+i — Xi < L for all i G Z, (59) 
and a number ttiq so that the ratio 

mo = ^+ 1"^" jg independent of ? G Z. (60) 

P-[Xi,Xi+l\ 

This assumption is not the weakest necessary assumption and we make it for clarity of exposition. 
The simplest example where this assumption is satisfied is periodic percolating flows for which (| 
or ( [5^ ) is satisfied. 

Then we have the following Theorem. 
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Figure 3: Streamlines ofu{x,y). 



Theorem 6 Let each of the sets be of the form = G [cj — hj, Cj + hj]}. Then under the 
assumptions made above, we have 



( 



{y)r > c 



+ 



1 



1 + mo 



-1 ^^+^ 



\ 



1 



-1 



J \u{p,9)\E,{p,e)^ 



(61) 



2 

h. 



1 + mg 



tE 1 + 



\u{p,e)\E^ip,e) 



dp 



^3 2 



for every 



T > tq = max 



K H + L 



,2 ' 



Here I = k/vq, L andniQ are as in (5i) and (6C), and the constant C in ( \6]\ ) depends on the function 
f{T) and the constants appearing in (54) and (St). 

Remark. Notice that the integrals on the right-hand side are independent of and give fluxes of the 
fluid through the middles of the tubes of streamlines. 

Proof. The proof of this theorem follows the steps of the proof of Theorem ^ for the shear flow. 
We will again utilize the differential inequality (^), and the expression ( p!7| ) for the bulk burning 
rate, as well as multiple averaging over regions bounded by the streamlines of the advecting velocity 
u{x,y). 
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Let us consider one region = {{p,6) : p G [cj — hj, cj + hj]}. Let us also denote the tube of 
streamlines = {{p,0) : p G [cj — 6,Cj + d],6 G {—00,00)}, 6 < hj and integrate (||) over the set 
Ds C Dp 

s 2 
J Ttdxdy - J dpu{p,0)Ei{p,0) + k J ATdxdy = ^J f{T)dxdy. (62) 



Ds -S 
We used here the relation 



Ds 





J u-VTdxdy = — J u{p,9)Ei{p,6,)dp 
Ds -S 



(63) 



which follows from the fact that is the tube of streamlines, and the boundary conditions 
Moreover, the quantity on the right side of (|63| ) is independent of 6 because u{x, y) is incompressible 
i). 

We first estimate the term that involves the Laplacian in (62). The following analog of Lemma 
^ holds: 
Lemma 6 

j d-i J d6 J AT{x,y)dxdy <C J [h]\VT\^ + f{T)]dxdy. 

hj /2-y Ds 

Proof. Notice that 

AT{x, y)dxdy = 



dO 



EidT 



E2dT, 



Ds 



dO 



dT dT 

u;{6,e)—i5,e)-coi-5,e)—{-6,i 



(64) 



with oj{d, p) defined in (^5|). Next, following the general procedure in the proof of Theorem ^ we fix 
9 £ M and average (|6^) over 5 G [hj/2 — 7, hj/2 + 7] with 7 G [0, hj/2], and then also average in 7. 
Then (pi) becomes 



hj /2 hj /2+7 



dj 



d6 [Tp{6, 0M5, e) - Tpi-5, 9)u;{-6, 9)] . 



(65) 



hj/2~-i 



We show how to estimate the first term in (p5|), with the second term treated in the same way. We 
integrate it by parts to get 



hj /2 hj /2+7 



d-i / d5Tf,{5,e)uj{5,e) 



, hi 



hj/2 



d'y 



?tj/2-7 

,h 



r(^ + 7, 0)0.(^+7,0) 



hj/2 /ij/2+7 

J d-/ J d6T{6,9)ujp{6,9) 

/ij/2-7 



hj/2 



d'J 



,hi 



,h 



n^+l,d) -l]u;U + 7,9)-{T{^-j,9)-l\ u;{^ - 7 



,hi 



,hi 



(66) 
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hjj2 hj/2+1 

J J d6{T{6,9)-l)u;p{6,e). 

/ij/2-7 



Consider the set of 6 such that 



We have for such 6 



hj/2 



hj/2 /ij /2+7 

j d-i j d5\Tp{5,e)u;{d,e)\<C j djV2^^j d5T^{d,( 



Next we look at 9 such that 



1/2 

< Ch] J Tjdp. (67) 



In this case for any pi, p2 G [0, /ij] we have 

\T{pue)-T{p2,e)\<^h, 



1/2 



Therefore, either r(/), 6*) > 1/4 for ah p, or |T(p,6')-l| > 1/4 for all p. Then we have, 1-T < C/(T), 
or T < C/(T), respectively. We use one of these bounds and the corresponding part of (|6^) to get 
for such 9 



hj/2 /ij/2+7 

J d-f J d5Tp{5,9)uj{6,9) 

hj/2-^ 

<C [ f{T{p,9))dp. 



hj/2 /ij/2+7 



<C J f{T{p,9))dp+- J dj J dSf{T{5,9)) 

hj/2-j 



(68) 



Now we put together the estimates (^^ and ( |68D to obtain for all 9 (z M: 

hjj2 ?tj/2+7 oo 

d-i f d6 [ d9Tp{6,9)uj{6,9) 



hj/2~^ -oo 



(69) 



< C 

< C 



h] j d9 j dpT^+ j d9 j dpf{T{p,9)) 



-oo 



-oo 



h]\VT\^ + f{T{x,y))\dxdy. □ 
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Similarly to the shear case, Lemma ^ and (|6^) imply the inequality 



< C 



H 



K / IVT 



H h] + P 



2 dxdy v"^ 



H 



dxdy 
H 



(70) 



where G{h,^) is defined as before by (27). An estimate similar to ( [70| ) holds in the regions D~ , 
where the flow is going backwards, except that the time derivative term in (^0|) enters now with the 
opposite sign: 



d: 



/ |^(|^|„(p,«)|£,(,,«) (71) 



< c 



n / ivn 



dxdy Vq 
H ^ 4^ 



f{T{x,y)) 



dxdy 
H 



D~ 



d: 



Let us choose the weights 



1 + mo 



1 + m, 



-1 ' 



(recall mo is defined by (pOj)) so that 

m+/i+[xi,a:i+i] = m_/i_[xi, 



(72) 



for any two points Xj, Xj+i of the partition (^) of the x-axis, similarly to what we did in the proof 
of Theorem Q (see (37)). In order to finish the proof of Theorem ^ we multiply equations (^) and 
(|7l| ) by m+ and m_, respectively, and add them. It remains now to estimate the time derivative 
term, and the following general Lemma provides us with the analog of Lemma ^ in the shear case. 

Lemma 7 Let be a rectangle Qq = Lx H , and let Qi^2 C rio be two open subsets o/r^o- Consider 
two continuous non-negative functions 4>i^2 ■ ^1,2 — ^ such that < 4>i^2ix,y) < C and 



dxdy(j)i{x,y) = J dxdy(t)2{x,y). 

Hi ^2 



(73) 



Let T : ^ M be a continuously differentiable function, < T < 1, then for any e > we have 



dxdy<j)i{x,y)T{x,y) - J dxdy(j)2ix,y)T{x,y) 
Q2 



(74) 



< C 



(L + H) 



e I dxdy\VT\^ + -J f{T{x,y))dxdy + J f{T{x,y))dxdy 

S7q Q,q 
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We postpone the proof of Lemma |^ till the end of this section. Using Lemma |^ in each rectangle 
Rxi,xi+i with e = ^, r^i = 1)+ n Rxi,x,+i, 0.2 = n Rx„xi+i, and functions ^1,2 given by 



d: 



we arrive at the analog of Lemma ^: 



dt 



^+ E ^Hih^+m ^^^P^ e)xn^ (P, Omip, e)dpdd (75) 



<^E 



{H + L)k f 2dxdy f (H + L)v()\ f j,,^, ^^dxdy 



i=l 



\VT{x,y,Ti 



H 



+ 1 + 



K 



f{T{x,y,Ti))- 



H 



n 



The rest of the proof of Theorem 



is completely analogous to the proof of Theorem ^. We average 
l"7|) and (|T^) to conclude the proof. □ 



in time according to (45) and use 

We now give the proof of Lemma |^. 
Proof. We define the measures ^1^2 by 

di^i,2{x,y) = (l)i^2{x,y)xn^,2^x,y)dxdy. 

Let A <Z he the set of points where T(x,y) > 5/8, and let the open set S C ^2 be such that 
I'liA) = i'2{Bi). Then we have 

di^i{x,y)T{x,y) - J du2ix,y)T{x,y) < C J dxdyf{T{x,y)) (76) 

+ j dui{x,y)T{x,y) - J dv2{x,y)T{x,y). 
A B 

Let us decompose further B = B' U B" , where B' = {{x,y) G B : T{x,y) > 3/4}. We also consider 
an open set A' C A such that vi^A') = i'2{B'), and write A = A' U A". Then we obtain 



J dnuiT — J dv2T = j dviT — j dv2T + j dviT — j dv2T 
A B A' B' A" B" 

and, moreover, Lemma ^ implies that 

dviT{x,y) - j dv2T{x,y) <C j dxdy f {T {x , y)) + J dxdyf{T{x,y)). 
A' B' A' B' 



(77) 



(78) 
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L 

Figure 4: Staircase for 1^ and J^. 

Therefore, we are done if ui{A!') = 1^2(13") = 0. Assume now that this is not the case. Then we may 
find a horizontal line h : y = yo and a vertical line I2 '■ x = xq such that 

|/i n A"\ > Y^iiA"), \l2 n B"\ > ^MB"), 

h rl 

where |S| denotes the one-dimensional Lebesgue measure. Moreover, we may choose subsets Qi C 
h n A" and Q2 <Z B" so that Qi = U^^^/fc and Q2 = U^LiJk are finite unions of intervals, and 
Cui(A") 

\Qi\ = IQ2I ^ • We may assume (possibly after subdividing into smaller intervals) that 

L -\- H 

N = M, and = \Jk\ for all k. Let us connect each pair of intervals Ik and Jk by perpendicular 
lines "staircase" as depicted on Figure ^. Notice that for every point {x,y) G Ik, T{x,y) > 3/4 while 
for every point {x' , y') G J^, T{x' , y') < 5/8. An argument directly analogous to that in the proof of 
Lemma shows that the following estimate holds: 



J f{T)dxdy J \VT\^dxdy > C\Ik 



Therefore we have for any e > 



^ j f{T)dxdy + e J \VT\^dxdy > ^^„,{A") 



> 



C 



L + H 



f2o 

dxdy(j)i{x,y)T{x,y) - J dxdy(p2ix,y)T{x,y) 
A" B" 



(79) 



25 



Equations (|7^[79|) show that 



dxdy(j)i{x,y)T{x,y) - J dxdy(j)2ix,y)T{x,y) 



(80) 



< C 



(L + H) 



e J dxdy\VT\^ + -J f{T{x,y))dxdy + J f{T{x,y))dxdy 



no 



The same proof shows that this bound holds for J^^ dxdy(j)2{x,y)T{x,y) — /^^ dxdy(j)i{x,y)T{x,y). 
This finishes the proof of Lemma 0. □ 

7 Examples of sub-linear growth of the bulk burning rate 

We give in this section examples of flows for which bulk burning rate grows sub-linearly in the 
amplitude of the advecting velocity. We do not try to identify the most general class of such flows 
but consider rather one simple family of flows of the form 



u{x,y) = UV^^ni{x,y) = ULy 



\ dy 



dx 



with the stream function 



"^mix, y) = cos""(7rx/L^) cos'^ {Try /Ly), m > 1 



(81) 



(82) 



periodic in x and y. The scalar U has the dimension of velocity. The structure of the level sets 
of the functions is clearly the same for all m. The period cell for these flows is the rectangle 
^ — [~^^^^] ^ [~"^>^^]> that consists of four smaller rectangles separated by separatrices 
^1 = 0. The normal component of u{x, y) is equal to zero at the boundary of the period cell of 
u{x, y), which slows down the burning as compared to percolating flows. This effect is quantified by 
the following Proposition. 

Proposition 1 Let T{x, y, t) be the solution of the reaction diffusion equation (Qj with either Neu- 
mann or periodic boundary conditions and fQj, respectively. Let u{x, y) be given by 
with Ly = H/2. Moreover, assume that the initial data TQ{x,y) has the property that TQ{x,y) = 1 
for X < xq, and TQ{x,y) = for x > xi. Then there exists a constant C > such that for U > vq 
we have 

{V)oo . , I \ / [7^2/(1+-) ^ 



(83) 



Proof. We will construct a function </>(x), independent of y (and hence satisfying both the Neumann 
and periodic boundary conditions (^) and (Q)), and L^^-periodic in x, such that the function <l>(x, t) = 
g-A(a-ct)0^^'j satisfies the inequality 



$t + u • V$ - kA$ > 0. 

4k 



(84) 



Moreover, the function (/>(x) will be positive, bounded, and bounded away from zero. Then maximum 
principle will imply that the solution T{x, y, t) of (|l|) with the Neumann or periodic boundary 
conditions satisfies the inequality 
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since it holds at t = for all A > because of our choice of the initial data. Then we will have 

(V^)oo < c (85) 

as in Lemma |^. Therefore our goal is to find a function (j){x) and A > so as to satisfy with as 
small c as possible. The function (t){x) should obey the inequality 

L(j) = — — 2k,(I)x — — ■ V4> + ui4> < 
A A 



with 



B = c — kX 



(87) 



We will define 4>{x) on the interval [—Lx/2,3Lx/2], and then extend it periodically to the whole 
real line. In order to make use of the fact that the x-component of u{x, y) is small near the lines 
X = —Lx/2,'iLx/2 we consider a smooth cut-off function x{x) defined as follows. Let rj{x) be a 
cut-off function 



r/(x) 



1, \x\ < 



2 \vo) 

and r/ decays monotonically from one to zero between those intervals, so that 



r/'|<C — , |r?"|<C — 



2a 



The exponent a > is to be chosen later. Then we define for x G [— 3Lx/2] 

1 



X{x) = ?? i— +o +^ 7 o ' 



Lx 



so that the two terms have non-overlapping support, and set 

(t>{x) = x{x) + (1 - X(a:))e^^ := x{x) + P{x) 

so that 

g-AL./2 <g3AL./2_ 

First we observe that since < CU (U/vq)""^'^ on the support of xix)i we have 
\Lx{x)\ = 

< C 



jX - 2kx - yX +uiX 



XLl \vq) 



XLx \vq 



Moreover, we have 
\L0\-- 



(1 - x)^^ - jx" + \uix: 



e^^' < C 



2a 



- + 



vo fU 



U 

+ t^o I — 



1— ma+« ' 



1— ma 



XLl \voJ ' XLx \vo 
We put together the bounds (^) and (p9|), and obtain 



e^^ + KA/3. (89) 



L(t> 



K fU 



XLl \vo 



2a 



+ 



+ VQ 



l~ma 



+ 



XLx \vo 



+ kX. (90) 
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Therefore the function (jj that we have constructed satisfies the inequahty (|86f) with the constant B 



given by the right side of (|90D. Using the definition (87) of and relation ( pq) we obtain then the 
following bound on the bulk burning rate: 



I 



2a 



I 



u 



1—ma 



+ 



1 



AL, 



4/A 



where / = k/vq is the laminar front width. Then we let A = l/L^ and a = 1/(1 + m), and get the 
estimate in Proposition 

One can see from the proof of Proposition [l| that it may be easily generalized to include cellular 
flows other than those of the form (^ll-|82|) . The relevant assumptions are similar geometric structure 
of the streamlines and the appropriate rate of decay of the normal velocity at the boundary of the 
period cell. 

The power ^ = 2/(1 + m) in (|8^ ) is probably not sharp, but Proposition |l] still shows several 
important points. Not all advection velocities with non-trivial ui lead to linear growth of the bulk 
burning rate in the advection amplitude. The presence of the closed streamlines appears to be crucial 
for sub-linear enhancement. Also, the exponent ^ may be made arbitrarily close to ^ = by taking 
m — > oo, and thus one can construct non-trivial flows for which the bulk burning rate grows slower 
than any given power of U/vq. 



A Homogenization regime 

Here we briefly present a simple direct application of the bound we obtained in Section |3|. We will 
consider a homogenization regime where the reaction is very weak, and investigate an effect of the 
periodic advection velocity in this limit. Let us consider the reaction-diffusion-ad vection equation 
(H) with the laminar velocity vq being small: vq — > jj'Vo, N ^ 1. The domain is then taken to be 
finite but very large: D^^ = ND, where D is some fixed region, such as a rectangle. Initial data 
varies on the large scale: 

{TN)t + n(x, t) ■ VTn = kATn + 7^/(?iv) (91) 

— — = on oDn 
on 

r^(x,o) = ro(-^), ^€Dn 



Then after rescaling x, t ^ A^x, N t the rescaled problem is 



,2 



{TN)t + Nu{N^, Nh) ■ VTn = kATn + ^/(Tn) (92) 

— — = on oD 

on 

T7v(x,o) =ro(x), ^eDN 

We assume that n(x, t) is periodic in x with period cell Q and in t with period r, and vanishes on 
the boundary of the cell C. Moreover, it is convenient to assume that D contains an integer number 
of cells, so that u{Nx, Nt) vanishes on the boundary dD. The bulk burning rate is given as before 

by 

D D 

The following Theorem may be established using the technique of [^. 
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Theorem 7 The family of solutions of {9i) converges strongly in L2([0,r] x D) to the solution 
T of the homogenized problem 



<,7^ + TzfiT) (93) 



dxidxj ' 4k 

df 

— — = on dD 
on 

T(x,0) =ro(x). 



The tensor is given by 



T 

^ J ds J dyui{y,s)9j{y,s) 



Q 

with 9i being the periodic solution of the cell problem 

dO- 

+ u(x, t) ■ VOi - KMi = -u,(x, t). 

at 

C 

Moreover, there exists a constant C such that ||Tjv — 7'||i2([o,r]x_D) ^ 

Since D is finite, Theorem ^ implies that the bulk burning rate V]y{t) V{t) = / Tt{s)ds. Let 

Jd 

us denote /c* the minimal eigenvalue of the symmetric part of the tensor k* , and set Vq = vq \/k^/K. 
This Theorem may be also applied to the front propagation problem in a finite rectangle with the 
boundary conditions (|8|), @. Arguments similar to Theorem |l| imply that the bulk burning rate 
for the homogenized problem obeys the lower bound 



no >C?;oW^(l-e-"^o*/2.) 
for times less than t^, ~ diamD ^ Therefore the bulk burning rate in the original unsealed variables is 

1 



increased from j^-vq to — yfc^/K, but remains of order 0{jj). The dependence of the tensor k* on 

the advection velocity u and diffusivity k is rather complicated. Some estimates for k* were obtained 
in and [^], and they may be applied to obtain the relevant bounds for V{t). This homogenization 
analysis is applicable only in the limit of very weak reaction at a fixed diffusivity. This is the case 
when the front width Iq = k/vq is much larger than the typical scale of variations of the turbulent 
velocity. 
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